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Abstract Rare event search experiments using liquid 
xenon as target and detection medium require ultra- 
low background levels to fully exploit their physics po¬ 
tential. Cosmic ray induced activation of the detector 
components and, even more importantly, of the xenon 
itself during production, transportation and storage at 
the Earth’s surface, might result in the production of 
radioactive isotopes with long half-lives, with a pos¬ 
sible impact on the expected background. We present 
the first dedicated study on the cosmogenic activation 
of xenon after 345 days of exposure to cosmic rays at 
the Jungfraujoch research station at 3470 m above sea 
level, complemented by a study of copper which has 
been activated simultaneously. We have directly ob¬ 
served the production of ^Be, ^°^Rh, ^^®Sb, and 
^^^Xe in xenon, out of which only ^^®Sb could poten¬ 
tially lead to background for a multi-ton scale dark mat¬ 
ter search. The production rates for five out of eight 
studied radioactive isotopes in copper are in agreement 
with the only existing dedicated activation measure¬ 
ment, while we observe lower rates for the remaining 
ones. The specific saturation activities for both samples 
are also compared to predictions obtained with com¬ 
monly used software packages, where we observe some 
underpredictions, especially for xenon activation. 

1 Introduction 

Liquid xenon (LXe) is used as detection medium in cur¬ 
rent and future rare event search experiments, such as 
direct dark matter [1,2, 3,4, 5] and neutrinoless double¬ 
beta decay searches [6,7,8]. It features a high scintil- 

^e-mail: lbaudis@physik.uzh.ch 
'^e-mail: alexander.kish@physik.uzh.ch 
'^e-mail: fpiastra@physik.uzh.ch 
e-mail: marc.schumann@lhep.unibe.ch 


lation and ionization yield [9,10], as well as a high 
radio-purity. Apart from long-lived double-beta emit¬ 
ters, such as ^^'^Xe, ^^®Xe, ^^^Xe and ^^®Xe, where only 
the decay of ^^®Xe has been observed so far, xenon has 
no unstable isotopes. However, the exposure to cos¬ 
mic rays during production, transportation and stor¬ 
age aboveground can produce instable radio-isotopes 
via nuclear activation processes. 

For a given exposure time and altitude, the acti¬ 
vation yield of materials can be predicted using soft¬ 
ware packages such as Activia [11] and Cosmo [12]. 
Calculations with both codes were performed for the 
XENONIOO experiment [1], however the predicted pro¬ 
duction rates were too high to be compatible with the 
measured background rates [13,14,15]. 

In this work we present the first dedicated exper¬ 
imental measurement of the cosmogenic activation of 
a natural xenon sample, which we exposed to cosmic 
rays for 345 days at an altitude of 3470 m. In order to 
provide a benchmark for our activation measurements 
and predictions, a sample of oxygen-free high thermal 
conductivity (OFHC) copper has been simultaneously 
exposed at the same location. Due to its very high pu¬ 
rity, and thus low radioactivity levels, OFHC copper is 
one of the most-frequently used materials to construct 
low-background detectors, and a well validated material 
in the software packages. 

The intrinsic radioactivity of the samples, initially 
stored underground at the Laboratori Nazionali del Gran 
Sasso (LNGS, Italy) for more than 1.5 y, was measured 
by means of a high-purity germanium q-spectrometer 
at LNGS, before and after the exposure to cosmic rays. 
Section 2 describes the samples and their handling, Sec¬ 
tion 3 the cosmic activation procedure and its mod¬ 
elling, and Section 4 provides details on the q-spectro- 
meter and the data analysis. The results of the mea¬ 
surement and the comparison with predictions are pre- 
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Table 1 Isotopic composition of natural xenon and copper [16]. The half-lives of the double-beta decay isotopes ^^^Xe, ^®'*Xe 
and are from the NuDat 2.6 database [17], the range for the double electron capture isotope ^^®Xe is a theoretical 

prediction from [18]. 


Xe mass number 
Abundance [%] 

Li /2 [y] 

124 

0.09(1) 

>1.6x1014 

126 

0.09(1) 
[5-12] X 1025 

128 

1.92(3) 

129 

26.44(24) 

130 

4.08(2) 

131 

21.18(3) 

132 

26.89(6) 

134 

10.44(10) 

>5.8x1022 

136 

8.87(16) 

2.165x1021 

Cu mass number 
Abundance [%] 

63 

69.17(3) 

65 

30.83(3) 









Table 2 Altitude, atmospheric depth, and vertical flux of protons and neutrons, as well as scaling factors to relate the flux to 
sea level values for the various locations in our study. 


Location 

Altitude [m] 

Atm. depth [g/cm2] 

Vert, flux [m 2s isr i] 

Scaling factor 

Sea level 

0 

1030 

2.6 

1.0 

Lauterbrunnen 

795 

954 

4.7 

1.8 

LNGS 

985 

936 

5.4 

2.1 

Jungfraujoch 

3470 

728 

29.0 

11.2 


sented in Section 5, while conclusions are drawn in Sec¬ 
tion 6 . 


2 Samples and preparation 

The xenon sample consisted of 2.04 kg research-grade 
xenon (impurities <10 ppm) from Praxair with nat¬ 
ural isotopic composition (see Table 1). It was con¬ 
tained in a 1 liter stainless steel bottle at a pressure 
of ^100 bar (at T ~ 20°C). To avoid the measurement 
results to be dominated by activation products in the 
bottle (to = 2.9 kg), two identical bottles were used for 
the 7 -screening and for the activation procedure. The 
first one always remained underground, where the gas 
transfer took place as well. Both bottles were evacuated 
and baked at >100°C for several days before being filled 
with xenon. 

The 10.35 kg copper sample consisted of 5 blocks of 
OFHC copper (purity 99.99%, isotopic composition in 
Table 1) from Norddeutsche Affinerie (now Aurubis). It 
came from the batch used to construct inner parts of 
the XENONIOO detector (sample 6 in [19]). Before each 
measurement, surface contaminations were removed in 
an ultrasonic bath filled with the acid detergent Elma 
Clean 60 diluted with deionized water. The sample was 
then rinsed and wiped with pure ethanol (> 98%), and 
stored under boil-off N 2 atmosphere for several days in 
order to let the ^^^Rn diffuse out and decay. 


3 Activation by cosmic rays 

The activation took place in a controlled manner at 
the High Altitude Research Station Jungfraujoch [20], 
at an altitude of 3470 m above sea level, from October 
31®*, 2012 to October 15***, 2013. The cumulated activa¬ 
tion time was 345.0 days. The initial transport from the 
underground laboratory of LNGS to the Jungfraujoch 
took less than 5 days and is neglected in the analy¬ 
sis. Following the activation, the samples were brought 
to Lauterbrunnen (795m), where they were stored for 
about 4 days before being transported by car to LNGS 
for the underground 7 -measurement. Transport to and 
storage at the LNGS aboveground laboratory (985 m) 
lasted about 1 day. The cool-down times, from the time 
when the samples were brought underground until the 
start of the measurement, were 2.5 days and 14.8 days 
for the xenon and copper samples, respectively. 

For comparison with theoretical predictions and with 
other measurements, we have to convert our high-alti- 
tude activation results into specific saturation activities 
at sea level. We thus calculate scaling factors that relate 
the nucleon cosmic ray flux at the various locations, as 
shown in Table 2. The atmospheric depth at different 
altitudes was obtained using the U.S. Standard Atmo¬ 
sphere 1976 model [21], shown in Figure 1(a), with a 
sea level value of 1030 g/cm^ [22,23]. The vertical flux 
of protons and neutrons was then calculated for a given 
atmospheric depth based on Ref. [24], see Figure 1(b). 
Due to their lower flux [25] and lower cross section [26, 
27] muons only play a sub-dominant role in activation 
processes and are neglected here. 
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Fig. 1 (a) Relation between atmospheric depth and altitude from the U.S. Standard Atmosphere Model [21]. We adopt a sea 
level value of 1030 g/cm^ [22,23]. (b) Vertical nucleon flux as a function of the atmospheric depth [24]. (c) The cosmic ray 
spectrum from Activia and Cosmo, based on parametrization from Refs. [30,31]. Only the total flux is used in the calculations 
as protons and neutrons interact similarly at these energies. 


The cosmogenic activation was predicted using Ac¬ 
tivia [11] and an updated and revised version of Cosmo 
[12,29]. The cosmic ray spectrum encoded in the pro¬ 
grams is based on the parametrization from Refs. [30, 
31] and is shown in Figure 1(c). It was sampled from 
10 MeV to 10 GeV in 10 MeV intervals: the lower bound¬ 
ary is chosen because the energy thresholds of the nu¬ 
clear excitation functions are above this energy. The 
cosmic ray flux above 10 GeV does not affect our results 
due to the exponential behavior of the energy spectrum. 
Both packages calculate the cross-sections of the rele¬ 
vant nuclear processes, such as spallation, fission, and 
evaporation, using semi-empirical formulae developed 
by Silberberg and Tsao [28], with identical parame¬ 
ters for protons and neutrons. To distinguish between 
spallation, fission, and nuclear breakdown reactions, in¬ 
cluding contributions from fast fragmentation processes 
and transition areas, these formulae are defined in sep¬ 
arate regions, depending on the masses of the target 
and product nuclides (light, intermediate and heavy iso¬ 
topes). Above a nucleon energy of ^3 GeV, the cross- 
sections are assumed to be energy-independent [11,28]. 

Our calculations of the production rates were per¬ 
formed for an exposure to cosmic rays at sea level. To 
convert into saturation activities, we assumed an acti¬ 
vation time of 100 y and no cool-down time. 

4 Measurement and data analysis 

Before and after activation, the intrinsic radioactivity 
of the xenon and copper samples was measured with 
the high-purity germanium y-spectrometer Gator [32], 
operated underground at LNGS. The detector has a 
background rate of 230counts/day in the 100-2700 keV 


interval. Due to the ~3500 m water equivalent of rock 
shielding from cosmic radiation, the hadronic compo¬ 
nent of the cosmic radiation is completely absent and 
the atmospheric muon flux is suppressed by 6 orders 
of magnitude with respect to the aboveground labora¬ 
tory [33]. Hence it is safe to assume that no significant 
cosmogenic activation took place once the samples were 
stored underground. 

To establish the detection sensitivity to the y-lines 
of the expected activation products as a function of 
measuring time, the xenon and the copper samples were 
measured pre-activation for 26.5 days and 34.3 days, 
respectively. The post-activation spectra were acquired 
for 11.5 days for xenon and for 4.0 days for copper. The 
specific activities of the produced radio-isotopes are in¬ 
ferred from the intensities of their most prominent full 
absorption peaks in the post-activation spectra. The 
information in the pre-activation measurements is not 
used in the analysis as no cosmogenic isotopes were 
present. Both, pre-activation and post-activation spec¬ 
tra are shown in Figures 2 and 3. 

A Bayesian method, described in Refs. [34,35,36], 
was employed to infer the activity of an isotope using 
all its y-lines with a sufficiently large branching ratio. 
The spectrum around each line is divided into three re¬ 
gions, their width is related to the energy-dependent 
resolution a of the detector which is given in [32]. The 
signal region is defined as ±3 (t around the mean po¬ 
sition of the full absorption peak. The background is 
inferred from the count rates in two control regions, 
-|-3tT above and —3cr below the signal region. The total 
count rate ys in the signal region is given by 

ys = m ■ A{e ■ BR) , (1) 

wl + wr 
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Fig. 2 Pre- and post-activation spectra of the 2.04 kg xenon sample. One can identify the ^^^Xe lines at 202.9 keV and 
375.0 keV, and the ^^®1 line at 388.6 keV. Other prominent lines are from radioactive contaminations in the stainless steel 
bottle containing the xenon (primordial 238^ 232 ^^^ chains, cosmogenic ®^Mn, ®^Co). 
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Fig. 3 Pre- and post-activation spectra of the 10.35 kg OFHC copper sample. The entire post-activation spectrum is dominated 
by the cosmogenic activation products. 


where m is the sample mass, A is the specific activity 
of the sample, , 8 l,r are the background rates in the left 
and right control regions, respectively, and ws.l.r are 
the widths of the three regions. Hence the background 
rate in the signal region is the interpolation of the back¬ 
ground rates in the control regions. The product of 
branching ratio BR and detection efficiency e of the 
7 -line is calculated in Monte Carlo simulations, based 
on a detailed implementation of the detector and sam¬ 
ple geometry in GEANT4 [37]. The decays of the radio¬ 
isotopes of interest are simulated using the G4Radioactive- 
Decay class, where branching ratios and directional cor¬ 
relations are taken into account. 

Considering N 7 -lines for the same isotope, the like¬ 
lihood of the model is 

N 

= n hsj) ■ /p(Cl, \I3lJ) ■ /p(C'r, \M, (2) 

k 

where the C's,,,Lfc,Rfc are the counts in the three re¬ 
gions for each line fc, and t is the measurement time. 
fp{C\fj.) = jCX e~^ is the Poisson distribution. 


We use flat priors for the parameters jk and 
and generate the posterior probability density function 
(PDF) of the specific activity A by Markov Chain Monte 
Carlo methods implemented in the Bayesian Analysis 
Toolkit (BAT) [38]. The marginalised posterior PDF is 
used to decide whether we can claim the detection of a 
line: if the global mode of the posterior PDF and the 
left edge of the shortest 68.3% credibility interval (C.I., 
green region in Figures 4(a)-4(d)) are positive, we cal¬ 
culate an activity, using the mode of the posterior PDF 
as its estimator and the shortest 68.3% C.I. as iltr un¬ 
certainty. If only the global mode is positive, but the 
left edge of the shortest 68.3% C.I. is zero, we report an 
upper limit as the signal is too weak to be determined. 
The upper limit is given as the 95.5% quantile of the 
posterior PDF (the right edge of the yellow region in 
Figures 4(e) and 4(f)). 


Because the half-lives of several examined isotopes 
are comparable to the integral time of the measurement. 
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Fig. 4 (a) Posterior PDF for the detected ^^^Xe activity, clearly visible in Figure 2. The green region represents the shortest 
68.3% C.I., which is taken as TIct uncertainty, (b) ^^®Sb, (c) ^^^^Rh and (d) ’^Be are detected as well, but with larger 
uncertainties, (e) For i 23 m'pg ^ signal is present, but an upper limit is reported since the intensity is too low to be determined, 
(f) An upper limit is reported for the activity of ®®Rh. The right edge of the yellow region represents the 95.5% upper limit 
(or the lower 95.5% C.I.). 


the measured mean specific activity A can significantly 
differ from the one at the beginning of the measurement 

Ao = A- (3) 


where r is the isotope’s mean lifetime, and t is the mea¬ 
surement time. 

The production rates at Jungfraujoch Aj are calcu¬ 
lated from the specific activities at the start of the mea¬ 
surement Ag, taking into account the activation time 
and the cool-down time For short-lived nuclides, we 
must also take into account the rather short time in¬ 
terval ft, where the samples were exposed to a reduced 
cosmic ray flux during storage and transportation at 
lower altitudes, leading to a lower activation rate At. 
The corrected specific activity Aj is 



ott/T 



(4) 


where r = AtjAj is the ratio of vertical nucleon fluxes 
at Lauterbrunnen and the Jungfraujoch. We combine 
the storage and transport time to ft = 5 days and as¬ 
sume an altitude of 795 m for the whole period. The 


specific activity of a sample of mass m after an activa¬ 
tion time ta is 

A{ta) = = P (l - . (5) 

m T V / 

The specific production (activation) rate P is equal to 
the specific saturation activity 

Kat = Mta -f oo) = P (6) 

for long activation times. In order to compare with 
other measurements and calculations, we convert A] 
to the specific saturation activity at sea level Agat by 
scaling it with the factor given in Table 2: 

Agat = Aj/11.2. (7) 


5 Results 

5.1 Copper 

The radioisotopes of interest were selected based on 
their half-life {A 1/2 >5 days), the expected production 
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Table 3 Results for the specific saturation activity Agat of natural copper at sea level, derived from our measurements of the 
cosmogenic activation. These are compared to our predictions from Activia and Cosmo, using semi-empirical formulae for the 
cross sections. We also compare to a measurement performed at LNGS [39] (scaled to sea level by a factor 2.1 and corrected for 
an evaluation error in case of ®'^Mn and ®®Fe [40]), to predictions with Activia [11] using (a) the same semi-empirical formulae 
and (b) the MENDL-2P database for the cross sections, to semi-analytical calculations [41] using cosmic ray spectra from (c) 
Ziegler and (d) Gordon et al, and to predictions using TALYS [42]. Deviations from our measured values beyond +lcr and 
— la are indicated by bold or italic font styles, respectively. 


Isotope 


Ti/2 


Copper: specific saturation activity at sea level Agat [z^Bq/kg] 



[days] 

This work 



Literature values 





Measurement 

Calculations 

Measurement 

Activia [11] 

Calc. [41] 

Calculation 




Activia 

Cosmo 

LNGS [39] 

a 

b 

c 

d 

TALYS [42] 

«Sc 

83.79 

27t“ 

36 

17 

25.2 ±8.6 

36 

36 

44 

31 

- 

> 

00 

15.97 

39t(® 

34 

36 

52± 19 

- 

- 

- 

- 

- 

54Mn 

312.12 

154t®f 

166 

156 

394 ± 39 

166 

145 

376 

321 

188 

59 Fe 

44.50 

^7+11 

49 

50 

57± 14 

49 

21 

75 

57 

- 

56Co 

77.24 

108+lt 

101 

81 

110 ± 14 

101 

163 

153 

231 

- 


271.74 

519+i“° 

376 

350 

860± 190 

376 

421 

1022 

858 

650 

5SCo 

70.86 

798t®8 

656 

632 

786 ± 43 

655 

441 

1840 

1430 

- 

®0Co 

1925.28 

340t«| 

304 

297 

1000 ± 90 

304 

112 

1130 

641 

537 

rate, and 
line at 

on their y-spectrum: we require 
> 60keV with BR>10%, which 

at least one 
is separated 

ment between measurement and predictions indicates 
the validity of our implementation of the cosmic ray 


from a prominent background line by more than 6 cr. 

Table 3 presents the results on the specific satura¬ 
tion activity at sea level Agat for copper, derived from 
our activation sample. The numbers are compared with 
our own predictions using the Activia and Cosmo codes, 
with another measurement performed at LNGS [39], 
and with additional predictions from the literature: these 
are based on the Activia package [11], semi-analytical 
calculations [41], and the TALYS 1.0 code [42,43]. De¬ 
viations from our measurement beyond the ±lcr level 
are indicated by bold (too high) or italic (too low) font 
styles. 

All cosmogenic radionuclides identified in the ac¬ 
tivated copper sample are well-known: they are pro¬ 
duced in spallation reactions from the stable isotopes 
®^Cu and ®^Cu [11,39,41], The overall agreement be¬ 
tween our measurement and the Activia/Cosmo calcu¬ 
lations is remarkable, with most of the predicted activ¬ 
ities within the 68% C.I. The highest specific satura¬ 
tion activity, about 0.8 mBq/kg, is measured for ®®Co 
which has a half-life of 71 d. This value is ~20% higher 
than the prediction. The only other isotope where we 
measure a higher saturation activity (~30%) than pre¬ 
dicted is ®^Co. For most of the isotopes the calcula¬ 
tions with Cosmo yield systematically lower activities 
(^10%), with the exception of (1.5x higher than 
Activia) and ®^Mn (2x higher). The general good agree- 


flux at different altitudes, and can be considered as a 
benchmark for the comparison of the measurement and 
predictions for the xenon sample. 

For five out of eight isotopes, our results for cop¬ 
per agree with the only other measurement available in 
the literature, which has been performed at LNGS [39]. 
For 54Mn, 57Co and 6°Co, we observe production rates 
which are (2.5 ± 0.6), (1.7 ± 0.5), and (2.9 ± 0.8) times 
lower, respectively. We obtain identical results with the 
Activia calculations in [11], when we use the same semi- 
empirical formulae to calculate the excitation functions 
(case a). The predictions using the MENDL-2P data¬ 
base [44] (case b) tend to underpredict the produc¬ 
tion rates. The semi-analytical study [41] predicts much 
higher production rates than observed for both tested 
cosmic ray spectra (cases c [45], d [46]). The TALYS- 
based work [42] yields reasonable values, which are 20- 
60% higher than measured, and show better agreement 
with our measurement than with the one of Ref. [39] . 


5.2 Xenon 

For the xenon sample, we present the measured satu¬ 
ration activities at sea level in Table 4, together with 
our Activia/Gosmo-based predictions. We also compare 
the results to measurement performed by the LUX Gol- 
laboration [47], and to predictions using the TALYS 
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Table 4 Results for the cosmogenic activation of natural xenon. The specific saturation activities at sea level are compared to 
our predictions based on Activia and Cosmo, to a measurement by LUX [47] and to a calculation using the TALYS code [42]. 
The half-lifes refer to the numbers used in the software packages. Deviations from our measured values beyond +1<j and —la 
are indicated by bold or italic font styles, respectively. 


Isotope 

Ti/2 

Xenon: specific saturation activity at sea level Asat 

[^Bq/kg] 

[days] 

This work 


Literature values 


Measurement 

Calculations 

Activia Cosmo 

Measurement 

LUX [47] 

Calculation 
TALYS [42] 


53.3 

37012" 

6.4 

6.4 

- 

- 

85Sr 

64.8 

< 34 

5.3 

4.6 

- 

- 

88 zr 

83.4 

< 52 

6.7 

4.6 

- 

- 

giniNb 

62.0 

< 1200 

5.6 

5.0 

- 

- 

®9Rh 

15.0 

< 120 

8.3 

8.2 

- 

- 

lOiRh 

1205.3 

14201®" 

16.6 

15.3 

- 

0.5 

llOm^g 

252.0 

< 49 

0.9 

0.8 

- 

- 


115.0 

< 55 

51 

47 

- 

- 

125Sb 

986.0 

590121® 

0.2 

13.5 

- 

0.5 

121mrpe 

154.0 

< 1200 

299 

276 

- 

135 

123m^g 

119.7 

< 610 

14.7 

14.4 

- 

140 

126j 

13.0 

1751^^ 

247 

247 

- 

- 

131J 

8.04 

< 190 

147 

170 

- 

- 


36.4 

187012" 

415 

555 

1530 ± 300 

- 

129mXe 

8.89 

< 8.7 X 10® 

238 

421 

1360 ± 250 

- 

131mXe 

11.77 

< 3.6 X 104 

251 

313 

1620 ± 370 

- 


5.25 

< 1.2 X 10® 

159 

196 

1140 ± 230 

- 

132Cs 

6.47 

< 120 

166 

164 

- 

- 


code [42]. To re-scale the LUX-numbers to saturation 
activities at sea level, we use the procedure described 
in Section 3 and the information provided in Ref. [47]: 
we assume that all xenon was activated at sea level for 
150 days, followed by 49 days (7 days) activation of 50% 
(50%) of the inventory at the SURF aboveground lab¬ 
oratory. The subsequent cool-down time underground 
was 90days (132days). The atmospheric depth at the 
SURF altitude of 1600 m is 881 g/cm^, corresponding to 
a vertical nucleon flux of 8.5 m“^s“^sr“^. This yields a 
conversion factor of 3.3, in agreement with the number 
given by LUX [47]. 

After activation, we have detected y-lines from the 
following isotopes: ^Be, ^°^Rh, ^^®Sb, ^^®1 and ^^^Xe. 
While our measurement and predictions for agree 
within the statistical and systematic errors, this is not 
the case for the light isotope ^Be, where the measure¬ 
ment is ^50 times higher than the prediction, and for 
^^®Sb, where we observe a (2900 ± 1200) times higher 
activity than predicted by Activia and a (44 ± 17) times 
higher activity than the Cosmo prediction. The ob¬ 
served production rate of ^^^Xe agrees with the mea¬ 
surement in the LUX detector [47] , however, the predic¬ 


tions are about a factor 4 too low. Our sensitivity did 
not allow us to detect the short-lived xenon isotopes 
i29“Xe, i^imxe, and ^^Sxe. This also holds for the var¬ 
ious other isotopes predicted by Activia and Cosmo, or 
by the study using TALYS [42]. We note that the cos¬ 
mogenic production rates for xenon isotopes predicted 
by Cosmo are systematically higher than the ones from 
Activia. 


6 Discussion and Conclusions 

We have carried out the first experiment to directly 
study the cosmogenic activation of xenon, which is em¬ 
ployed as target and detection medium for rare event 
searches, and the second dedicated measurement on the 
activation of OFHC copper, often used as ultra-pure 
detector construction material. Both samples were ac¬ 
tivated by cosmic rays in a controlled manner at the 
Jungfraujoch research station (3470 m) to maximize 
the incident nucleon flux. The activation products were 
measured with a low-background germanium detector 
in the Gran Sasso underground laboratory. We have 
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compared the measurement results to our own predic¬ 
tions obtained with the Activia [11] and Cosmo [12] 
software packages, as well as to other measurements 
and predictions found in the literature. 

From the eight detected isotopes in the OFHC cop¬ 
per sample, the production rates of five agree with a 
measurement performed at LNGS in 2009 [39], which is 
the only dedicated measurement available in the litera¬ 
ture. We measured up to a factor of ~3 times lower val¬ 
ues for the other three isotopes, see Table 3. While the 
agreement with our Activia/Cosmo predictions is satis¬ 
factory, supporting the validity of our comic ray model 
for the xenon study, discrepancies with other published 
calculations are present. 

After cosmic activation of the xenon sample, we de¬ 
tected hve isotopes, four of which were not directly mea¬ 
sured before: ^Be, ^°^Rh, ^^®Sb, The measured 
saturation activity of ^^^Xe agrees with the measure¬ 
ment by the LUX collaboration [47]. In general, the 
agreement with predictions from Activia and Cosmos 
is unsatisfactory, the production rate for only one iso¬ 
tope, is calculated correctly. The predicted rates 
for most isotopes are too low by a factor of a few to 
^10^, which is potentially worrisome, as some of these 
have rather long half-lives (3years). This conclusion 
also holds for the short-lived xenon isotopes ^^®™Xe, 
isim^g and ^^^Xe which we could not detect in our 
study: all predicted numbers are significantly lower than 
the values measured by LUX [47]. 

The main sources of background in multi-ton scale 
LXe-based dark matter detectors are the ones which are 
uniformly distributed in the target, such as ^^^Rn and 
®®Kr, or interactions from solar neutrinos [50]. Among 
the observed cosmogenic activation products of xenon, 
only ^°^Rh and ^^®Sb have half-lives which are long 
enough to affect an experiment with a foreseen op¬ 
eration time-scale of ^5 years. Both elements have a 
relatively high electronegativity and could possibly be 
removed by the xenon purification systems that use 
hot zirconium getters. However, this has not yet been 
demonstrated experimentally. The decay of ^°^Rh does 
neither produce low-energetic electrons nor low-energetic 
7 -rays not accompanied by a prompt (<lns) second 
signal, which shifts the initial low energy signal to high 
energies by pile-up. We therefore conclude that ^°^Rh 
will not lead to dangerous background. The situation 
is different for ^^^Sb, where 13.6% of the beta-decays 
end up in a long-lived excited state of the ^^®Te daugh¬ 
ter (ri /2 = 57.4 d) and are therefore unaccompanied by 
a subsequent y-ray. The low-energy tail of these elec¬ 
trons will lead to single-scatter electronic recoil back¬ 
ground, while all other decay paths do not contribute, in 
agreement with the Geant4 prediction. Assuming that 


Sb is not removed by the purihcation system, and us¬ 
ing the activation times quoted by LUX [47] together 
with our measured activation rate, a background con¬ 
tribution of (3.0 ± 1.3) • 10“^ events • keV“^ • kg“^ • d“^ 
would be expected. This rate is too large compared to 
the published total background level of (3.6 ± 0.4) • 
10“^ events • keV“^ • kg“^ • d“^ [47]. We therefore con¬ 
clude that the true activation rate is either close to the 
lower end of our quoted credibility interval, or Sb is 
removed by the getter or plates out at surfaces. 

In addition, the chemically inert noble gas isotopes 
could also affect next-generation dark matter searches. 
We have only observed ^^"^Xe, which has a relatively 
short half-life (ri /2 = 36.3 d), similarly to other cosmo- 
genically produced Xe-isotopes found in [47] . These iso¬ 
topes will not pose a problem for next-generation LXe 
detectors, as their contamination will be reduced by 
an order of magnitude just after a few months of stor¬ 
age below ground. However, we note that for multi-ton 
scale LXe detectors such as LZ [48], XENONnT [49] and 
eventually DARWIN [5,50], the cosmogenic production 
of radioactive isotopes by the muon flux present at the 
underground location must be studied as well, as in this 
case also the decays of the short-lived isotopes might 
provide a source of backgrounds. 
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